VOLUME  I 

PERFORMANCE  FLIGHT  TESTING 


APPENDIX  F 
DERIVATIONS 


F.l  DEElIVATiasr  OF  EULER  EQUATICa^  (Chapter  2) 


Starting  with  Newton's  Second  Law 

F  =  ma 

for  an  eleanent  moving  along  a  streamline  with  velocity  V  in  the  X  direction 


FIGURE  F.1.1.  FORCES  ON  A  FLUID  PARTICLE 

EF  =  pdydz  -  (p  +  ^  .dx)  dydz 
or 

^  (dxdydz) 

mass  =  (volume) (density) 

i.e. , 

m  =  p  (dxdydz) 

^  dt  dx  dt  dx 
.*•  -  ^  (dxdydz)  =  p  (dxdydz)  V  ^ 


i 


F.2  DERIVATIOJ  OF  STANDARD  ATMOSPHERE  RELATIONS 


(Chapter  5) 


Beginning  with  Equation  5.4  fron  5.7 


G  ^ 
R  T 


Since  T  = 
a 


Therefore, 


f(H)  below  36,089  ft  due  to  the  lapse  rate,  L,  we  can  substitute 


Equate  F.2.1  and  F.2. 2,  the  two  sides  of  Equation  5.4 


In 


(■ 


H 


) 


F.2 


DERIVATICN  F.2  (Continued) 


or 


=  (1  -  K^H)^2 


Finally,  since 


“  ^a'^a 


we  can  get  the  expression  for  p  from 


(1  -  Kj^H)^2 

"  (1  -  K^H) 


K_-l 

a  =  (1  -  K^H) 


G/RL 


Similar  derivations  can  be  done  for  H  >  36,089  ft 


( 


F.3  DERIVATia^  OF  SPEED  OF  SOUND  (Chapter  6) 


FIGURE  F.3.1.  PRDPAGATKU  OF  A  SMALL  PRESSURE 
PULSE  IN  A  FRICnONLESS  PIPE 


With  the  piston  stationary  the  velocity  in  the  pipe  is  zero,  and  the 
fliiid  has  seme  density,  pressure,  and  terperature.  Displacing  the  piston  to 
the  dotted  position  (an  infinitely  small  distance)  causes  a  small  pressure 
pulse  (dP) .  This  small  dP  travels  down  the  pipe  at  a  speed  defined  as  the 
speed  of  somd,  a,  for  the  fluid  in  the  pipe.  This  speed  of  an  infinitely 
small  pressure  pulse  is  also  known  as  acoustic  velocity.  Behind  the  pressure 
piiLse  the  air  has  a  small  dV  because  the  piston  has  displaced  the  fluid. 
While  not  obvious,  this  fact  can  be  verified  experimentally.  To  analyze  this 
situation  mathematically,  a  coordinate  change  will  be  made.  Such  a  change  is 
ccinmon  in  almost  all  fluid  mechanics  and  aerodynamics  courses  and  is  the  basis 
for  wind  tunnel  testing.  The  math  model  in  Figure  F.3. 2  transforms  to  that 
shown  in  Figure  F.3. 3.  Notice  that  the  pressure  wave  in  Figure  F.3. 3  is 
stationary. 


F.4 


DERIVATION  F.3  (Continued) 


T  +  dT 


FIGURE  F.3. 2.  MOVING  WAVE  FIGURE  F.3. 3.  STATIONARY  WAVE 

Evaluating  the  monentum  equation  across  the  pressure  pulse  frcm  to  in 
Figure  F.3.3 

dP  +  pVdV  =  0 

dP  =  P2“^l  ^  P  +  dP-P  =  dP 

dV  =  "^2  ~  '^1  ^  a-dV-a  =  -dV 

Substituting 

dP  +  pa  (-dV)  =  0 

dP  =  padV  (F.3.1) 

Evaluating  the  continuity  equation  across  the  pressure  pulse  from  to 

gives 


pAV  =  Constant 
(pAV)2  =  (pAV)^ 

(p  +  dp)  (a  -  dV)  =  pa 


pa  -  pdV  +  adp  -  dpdV 
dp  =  (p/a)  dV 


pa 


(F.3. 2) 


DEEtCVATICN  F.3  (continued) 


Dividing  Equation  F.3.1  by  F.3. 2 


dP  _  padV 
dp  ~  (p/a)  dV 


F.6 


F.4  DERIVATICN  CF  STEADY  FLOW  ENERGY  EQUATION  (Chapter  6) 


FIGURE  F.4.1.  CGNTRDL  VDIUME  FOR  ENERGY  BALANCE 

de  =  dq  -  dW  (First  Law  of  Ihennoc^amics)  (F.4.1) 

dW  =  Work 


dWj,  E  Shaft  work 
J\ 

dq  H  Heat  transfer 


de  H  Energy  (Capacity  of  flow  to  do  work) .  Several  types 
of  energy  are  possible. 

First  kinetic  energy: 

KE  =  i  mV^ 
d(KE)  =  mVdV 


or  for  a  unit  mass 


d(KE)  =  VdV 


F.7 
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DERIVATIOT  F.4  (Continued) 

A  secxjnd  type  of  energy  is  potential  energy: 


d(PE)  =  gdz 

Internal  energy  is  energy  due  to  random  motion  of  the  molecules  of  the  fluid. 
Cjonsider  the  two  fluid  molecules  in  the  acccarpanying  sketch. 


O 

1 

Vi=V, 


FIGURE  F.4. 2.  MOLECULAR  MOTICN 

The  kinetic  energy  of  each  molecule  is  the  same,  because  the  average 
velocities  and  V2  are  equal.  However,  molecule  2  obvioxjsly  has  more  randan 
energy  and  therefore,  more  energy.  This  type  of  energy,  called  internal 
energy,  u,  is  solely  dependent  on  absolute  temperature .  This  type  of  energy 
is  a  fluid  property. 

Expansion  and  flow  work  are  separated  fron  shaft  work  and  will  be  written 
on  the  energy  side  of  Equation  F.4.1  by  convention.  E3q>ansion  work  is  the 
work  done  in  order  to  expand  or  corpress  a  unit  mass  of  matter 

dCEjpansion  Work)  =  Pdv 

Flow  work  is  the  work  done  in  moving  the  unit  mass  of  matter 
D(Flow  work)  =  vdP 
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DERIVATIC8SI  F.4  (Continued) 

Now,  substituting  all  of  these  energy  and  work  definitions  into  Equation  F.4.1 
VdV  +  gdz  +  du  +  Pdv  +  vdP  =  dq  -  dW^  (F.4.2) 

gdz  +  VdV  +  du  +  Pdv  +  vdP  =  dq  -  dW^ 


But 


Pdv  +  VdP  =  d  (Pv)  =  d  || 


VdV  +  du  +  d(^)=  gdz  =  dq  -  dW^ 


K 


(F.4.3) 


Itecalling  the  definition  of  enthalpy 


h  =  u  V  I 


dh  =  du  +  d 


(I) 


Substituting  into  Equation  F.4.3 

gdz  +  VdV  =  dh  =  dq  -  dW^^ 


Integrating 


g  -  =  h2  -  hj  + 


V^  -  V^ 

-V^tg(Z2 


^l) 


(F.4.4) 


(F.4. 5) 


Potential  energy  change  in  our  problems  can  be  considered  negligible  vhen 
compared  to  the  kinetic  and  internal  energy  changes  of  our  system.  In  the 
case  of  an  adiabatic  process  with  no  shaft  work,  Equation  F.4.5  reduces  to 


V?  -  V? 


or 


v^ 


h  +  —  =  constant 


(F.4. 6) 


F.9 


(Chapter  6) 


F.5  RELATICNSHIP  BETWEEN  M  AND  M 


Starting  with  Equation  F.4.6  in  Derivation  F.4 


V 

h  +  =  constant 


Substituting  for  h 


2 


+ 


CpT  = 


tCpT* 


Using  =  yRT  a*^  = 


Siibstituting  into  Equation  (F.4.6) 


For  local  sonic  conditions  M  =  1.0  =  —  /.  V*  =  a* 

a* 

a^  _  a*^  a*^  _  (y  -  D  a*^  +  2a*^ 

2  Y-1  2  y-1  2(y-1) 


F.IO 


DERIVATIOSF  F.5  (Continued) 


Dividing  by 


2 


+ 


V2 


a*^  (y  -  1  +  2)  ^  a*^  (y  +  1) 
2{y  -  1)  2(y  -  1) 


But 


So 


1  I  1  -  1  (y  +  1) 

2  m2  (y-1)“m*2  2  (y-1) 


(F. 


Solving  for  M'^ 


Or 


1  ^  Y  +  1  _  1  ^  Y  +  1  -  M*  (y  -  1) 

m2  (y-1)  M*2  2  (y  -  1)  2  2  (y  -  1)  M*2 


m2=-L^ 


2  2 
M* 


[l- 


Now  solve  F.5.1  for  M*"' 


M*2  =  (V  +  1)  ^ 

m2  (y  -  1)  +  2 


F.ll 


DERIVATIOT  F.6  (Continued) 

From  continuity  AV^  =  (^-^2  *^1^1  ~  ^2^2 


(F.6. 4) 


From  perfect  gas  law  P  =  pRT  or  ■=  =  p  Substitute  into  F.6. 4 


P  P 

1  2 
V  =  —  v„ 

1  RT2  2 


T  P  V 

2  2  2 


M  =  -  ;  V  =  Ma 

d. 


T2  P2M2 


1  YYPTi 


h  ^ 

P  M 


■^/yRT 


(F.6.5) 


Equate  (F.6. 3)  and  (F.6.5) 


(F.6. 6) 


Use  mcmentum  Equation  dP  +  PVdV  =  0  and  frcsn  continuity  of  pV  =  constant 

P  +  PV^  =  const 


F.13 


DERIVATION  F.6  (Continued) 


^1  +  pf 

^2  ^2 

V  =  m-^yRT 

P 

=  M 

—  D  4- 

1  YRT2  ^2 

^  ^  RT^ 

“  Jr^  *r  JYL 

2  RT2 

Factoring  each  side 

II 

1 - 1 

T 

+ 

_IU 

r— 1 

.  p,  [1  + 

^“fj 

^2 

1  +  ym^ 

■ 

1  +  YM2 

Substitute  into  Equation  F.6. 6  and  rearrange 

1  + 

1  + 

2 

Solving  for  M2 


Ml  Vl 


"2 


2  2 
M  +  — = — 
1  Y  -  1 

-  1 

Y  -  1  1 


(F.6. 6a) 


(F.6. 7) 


(F.6. 8) 


Substituting  Equation  F.6. 8  into  Equation  F.6. 6a  and  Equation  F.6. 3  gives 
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DERIVATIOSI  F.6  (Continued) 


.  .gjr......  ^  ~  ?■ 

Y  +  1  1  Y  +  1 


1  -  Y  +  2y 
__ 


^2 

2y  .,2  Y  -  1 

Y-  1  . _2_ 

'^1  “ 

Y  +  l‘*l  Y  +  1 

J  ^  (y  +  1)  (1^)_ 

P2/P1  can  be  similarly  attained  with  lots  of  algebra. 

P2  (y  +  1) 

•"l  2  +  (y  -  1) 


(F.6. 9) 


(F.6. 10) 
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F.7  shocks  in  SUPEE^SCMSIIC  FLOW;  (Chapter  6) 

A  shock  has  been  described  as  a  discontinuity  between  supersonic  and  sub¬ 
sonic  flow.  Nothing  has  been  said  concerning  the  conditions  under  vdiich  it 
can  or  cannot  occur.  First,  the  existence  of  a  shock  wave  must  be  physically 
jtistified,  and  then  the  conditions  that  must  exist  before  a  shock  will  form 
must  be  determined. 

It  will  be  stown  that  a  shock  is  a.  discontinuity  between  supersonic  and 
subsonic  flow  and  it  will  be  shown  that  a  shock  can  OJLY  occur  viien  flow  goes 
frcm  supersonic  to  subsonic  conditions. 


NORMAL SHOCK 

FIGUEE  F.7.1.  NORMAL  SHOCK 

The  mass  flow  rate  through  a  shock  is  a  constant,  i.e.,  no  mass  is  added 
or  destroyed  by  the  shock,  and  the  cross-sectional  area  through  the  shock  is 
assumed  constant.  Therefore,  the  continuity  and  mcmentum  equations  may  be 
written 

m  =  pVA  =  constant  (F.7.1) 

m  =  ~  constant 


dP  +  pVdV 


0 


(F.7. 2) 


DERIVATION  F.7  (Continued) 

Evaluating  Equation  F.7. 2  at  Stations  1  and  2  (Figure  F.7.1)  and  substituting 
m  for  pV  gives 


^2  -  Pi  =  n,  Vi  -  V2 


Dividing  the  monentum  equation  by  the  continuity  equation: 


P  P 

^2  1 


P2  ^2  Pi 


Vi  -V2 


(F.7. 3) 


2  yP 

Multiplying  Equation  F.7. 3  by  y  and  substituting  a  =  -^ 


=2 

a2 

^2 

fl 

(F.7. 4) 

^2 

^1 

'1  2' 


Writing  the  energy  equation  for  a  point  in  the  free  stream  and  at  local  sonic 
conditions ; 


CpTT  =  CpT  + 


1/2 


c  T*  +  1/2  V* 
P 


(F.7. 5) 


Substitute  the  following  values  for  c^T  and  V*  into  Equation  F.7. 5: 


c„T  =  c„ 
P  P 


(?) 


R 


c  -  c 
P 


T  = 


Y  -  1 


c  T*  =  ---*Y 

P  Y  -  1 


V*  =  a*,  since  M  =  1 


F.17 


DERIVATION  F.7  (Continued) 


and  solving  for 


Y  +  1 
2 


(F.7. 6) 


Substituting  this  equation  into  Equation  F.7. 4  and  rearranging  (since 
a*2  =  a*j^  =  a*) 

*2 

— -  (V^  -  V2)  =  ^1-^2  (F.7. 7) 


If  equals  Equation  F.7. 7  has  a  trivial  solution,  i.e. ,  0  =  0,  or 

that  no  velocity  discontinuity  exists  in  the  flew.  It  is  an  experimental  fact 

that  V,  does  not  equal  V.,  across  a  shock  and  that  shock  waves  are  present 
^  ^ 

under  certain  flow  conditions.  Dividing  by  (V^  -  V^)  admits  that  there  is  a 
velocity  discontinuity  in  the  mathematical  flow  description  and 


a* 


2 


V2 


(F.7. 8) 


a*  is  the  speed  of  sound  at  local  sonic  conditions  and  can  be  shewn  to  be 
a  constant  throu^  the  shock.  The  shock  is  assumed  to  be  an  adiabatic 
process,  therefore 


and 


T* 


constant 


(F.7. 9) 


DERIVATION  F.7  (Continued) 


Therefore 


T*i  =  T*2 


=  a*2 


From  Equation  F.7. 8  it  can  be  seen  that  if  is  greater  than  a*,  then 
must  be  less  than  a*  in  order  for  the  equality  to  hold.  This  can  be  written 


cind  from  definition  of  M* 


M*^  = 


(7.10) 


Equation  F.7. 10  shows  that  if  is  greater  than  1.0,  then  must  be 
less  than  1.0,  i.e. ,  if  =  2.0,  then  M*2  =  .5. 

If  there  is  a  velocity  discontinuity  in  the  flw,  then  the  velocity  on 
one  side  of  the  discontinuity  must  be  subsonic  and  on  the  other  side  must  be 
supersonic.  This  relationship  between  and  M*2  gives  no  insight  as  to 
which  side  of  the  shock  is  subsonic  and  which  side  is  supersonic. 

Supersonic  side  of  shock: 

Next  it  must  be  established  vhich  side  of  the  shock  must  be  supersonic. 
Experiments  have  proven  that  a  shock  occurs  only  vhen  the  upstream  Mach  is 
greater  than  1.0,  but  vhy?  In  answer  to  this  question,  an  equation  for  the 
change  in  entropy  has  been  written 


dT  „  dP 
dS  =  Cp  —  =  R  — 


F.19 


DERIVATICN  F.7  (Continued) 

Integrating  and  rearranging  this  expression 


R 


c 

_E 

R 


In  T 


-In  P 


(F.7. 11) 


where  1  and  2  refer  to  the  stations  upstream  and  downstream  of  the  shock  wave. 
Evaluating  this  equation  with  the  stagnation  properties  at  Stations  1  and  2. 


AS  _  _P  2  ,  ^2 

-  =  V 

1  1 


and  since  In  1  =  0 

2  1 


R 


(F.7. 12) 


The  second  law  of  thermocfynamics  states  that  entropy  may  only  increase. 

Therefore  AS  in  Equation  F.7. 12  is  positive. 

In  P  /P  must  be  negative  or  P_  /P_  <1.  If  an  equation  relating 

2  1  2  1 

P„  /  P„  to  M,  can  be  derived,  the  question  of  vdiich  side  of  a  shock  wave  is 
2  1  ^ 

supersonic  will  be  solved. 

Using  Equation  F.7. 11  to  equate  the  entropy  change  in  terms  of  free  stream 
conditions  to  the  entropy  change  in  terms  of  stagnation  conditions; 


-  In 
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DERIVATICM  F.7  (Continued) 


F.21 


Notice  that  vAien  M  =  1,  AS  =  0.  This  is  the  case  of  the  isentropic 

somd  v»ave  or  veakest  possible  normal  shock  since  changes  across  it  are  so 
small  that  no  entropy  change  is  produced. 

F.8  LINEAR  THIN  WING  TEIEORY  (ACKERET  THEORY)  (Chapter  6) 

To  develop  the  Ackeret  Theory,  the  following  must  be  satisfied; 

1.  Geanetric  and  trigonetric  flow  relations 

2.  Conservation  of  mass 

3.  Conservation  of  mcmentum 

Geanetric  and  Trigonetric  Relations; 

The  geanetry  of  an  expansion  flow  is  shown  in  Figure  F.8.1.  The 
flow  for  a  conopression  Mach  wave  is  exactly  the  same  except  d6,  dV,  and 
dVjj  are  negative.  Thus  the  expansion  case  equations  will  be  valid  for 
the  carpression  case  if  the  signs  are  reversed.  Only  the  expansion  case 
equations  will  be  developed. 


FIGURE  F.8.1.  EXPANSION  FLOW 
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DERIVATION  F.8  (Continued) 

Since  d6  is  small  then 

dV  =  change  in  magnitude  of  V 


N  Sin  p 


dV^  sin  p 


also 


d6  = 


dVj^  cos  p 
__ 


Substituting  Equation  F.8.1  into  Equation  F.8.2 


d6  =  —  ^ 

V  sin  p 


dV  1 


V  tan  p 


FIGURE  F.8.2. 


Frcm  definition  of  sine  and 


sin  p  = 


1 

M 


Figure  F.8.2  can  be  constructed.  Frcan  Figure  F.8.2 

1 


tan  p  = 


7 


-  1 


(F.8.1) 


(F.8.2) 


(F.8.3) 


(F.8.4) 


F.23 


DERIVATICN  F.8  (Continued 


and 


or 


d6 


V 


V 


(F.8.5) 


Conservation  of  Mass: 

If  we  consider  a  constant  area  across  the  Mach  wave  then  fran  the 
conservation  of  mass  equation 


^IN  ^2  ^2N 


and 


(F.8. 6) 


or 


P2  =  Pi  +  dp 


^2N  = 


Pi  -  (p^  +  dp)  (V^  +  dVj^) 


(F.8. 7) 


Siiiplifying  Equation  8.7  and  dropping  dpdV^  as  insignificant 

dp  +  p^  dV^  =  0  (F.8. 8) 

Conservation  of  Mcmentim: 

Parallel  to  the  Mach  waves  there  are  no  pressure  differential  forces 
and  thus  no  manentum  flux  change  parallel  to  the  Mach  waves  and  the 
conservation  of  monentim  equation  beccmes 

^^1  ^IN^  ^IT  ^  ^^2  ^2N^  ^2T 


(F.8. 9) 


DERIVATIOJ  F.8  (Continued) 


but  from  conservation  of  mass 


'^1  ^2  ^2N 


(F.8. 6) 


^IT  =  V2T 

There  is  a  pressure-  differential  normal  to  the  wave  and  by  Newton's  second 
law,  this  pressure  differential  per  unit  area  must  equal  the  rate  of  change  of 
mcmentum 


<^1  ""in)  ""in  -  <^2  V  ''2N  =  ^2  -  '’l  <^-8.10) 


^IN  ^2  ^2N 


(F.8. 6) 


^2N  =  ^IN-^'^^N 


(F.8. 11) 


^2  =  +  dP 


(F.8. 12) 


Substituting  Equations  F.8. 6,  F.8. 11,  and  F.8.12  into  Equation  F.8. 10  gives 


V,.,  dV-,  +  —  =  0 

IN  N 


<3V„  = 


“l  ''iH 


(F.8. 13) 
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DERIVATION  F.8  (Continued) 


Now  we  will  conbine  the  relationships  of  the  three  previous  sections. 
Si±»stituting  Equation  F.8.13  into  Equation  F.8. 8  gives 


■  “l  ^ 


dP 


V. 


=  0 


IN 


rearranging  gives 


dp 


V,  =  ^  =  a' 
IN 


but 


a2  = 


yRT^  =  Y 


Thus 


V. 


2  _  ][!l 

IN  p. 


or 


yp 


Pi  = 


V. 


1 

2 

IN 


(F.8. 14) 


Substituting  Equation  F.8.14  into  Equation  F.8.13  and  rearranging  gives 


YPn 


dP  = 


V. 


IN 


yPi 


(F.8. 15) 


substituting  Equations  F.8.1  and  F.8. 3  into  Equation  F.8.15  gives 


yP. 


dP  =  - 


V. 


M  dV 


IN 


Multiply  right  side  by  V/V  and  substituting  a  for  gives 


dP  =  -  yP^  ^ 
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F.9  DERIVATICN  OF  THRUST  EQUATION  FOR  A  TURBOJET  (Chapter  7) 


The  following  equations  from  thermodynamics  and  physics  will  be  required 
and  are  presented  here  for  quick  reference. 

Frcm  gas  (^mamics  for  one-dimensional  gas  flew 


TOTAL 

^STATIC 


1  + 

y 


(F.9.1) 


'TOTAL 

^STATIC 


(F.9. 2) 


For  an  isentropic  process 


(F.9. 3) 


Mach 


General  thrust  equation 


(F.9. 4) 


5  <^10 -V 


(F.9. 5) 
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DERIVATICN  F.9  (Continued) 


^20 

—  ’  '••II 

T32 

1 

T43 

T54 

^105 

TT  20 

TT  32 

TT  43 

•"■54 

105 

1  : 

1  : 

i  4 

i  i 

>  9 

FIGURE  F.9.1.  STATION  DESIGNATIONS  AND  PARAMETER  DEFINITIOJS 


Define 


"^20  ■  K 


'^32  T„ 


^"32  P„ 


'^43  T„ 


""43  P„ 


T  =  — 
54  T„ 


54  P„ 


105  T„ 


105  P„ 


NOTE: 


'^105  ^^43  ^105  “  ^ 
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DERIVATICN  F.9  (Continued) 


Required:  ~  ^  (tertperatures) 


Solution: 


i^ly  Equation  F.9.1 


^10 


m  rp  fp  fn  <T>  fp 

^10  ^10  32  Us  U4  U05 


and  note  =  1 


Apply  Equation  F.9. 6 


^TIO  ^10 


1 


-  ^o  ^20  ^^32  Us  ^^54  ^105 


and  note  ~  ’’^105  ~  ^ 


Assume  the  nozzle  e}q)ands  gas  to  ambient  pressure  so  P,  „  =  P 


Equation  F.9. 7  yields 


1 


-  ir2Q  ^54  y 


(F.9. 8) 


From  Equation  F.9. 6 


m  m  m  0-1  m  m 

^10  ^  20  20  32  43  54 

To  ■  (l 


(F.9. 9) 


Combining  Equations  F.9. 8  and  F.9.10 


^^20  ^^32  ^43  ^^54 

h  P  P  I=i' 

UO  U2  ^54  Y 


(F.9.10) 
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DERIVATION  F.9.  (Continued) 


From  Equation  F.9.3 


20  Y 


=  T  If 

^20  '  32  Y 


irk  =  t:  ^  III  ^ 

^32  '  ^^54  Y 


=  T 


54 


Sxibstituting  these  e3<pressions  into  Equation  F.9.10 

T, 


10 


=  T 


43 


Dividing  Equation  F.9.9  by  Equation  F.9.11 


“  "^20  "^32  '^54 


Solve  for  M, 


10* 


M, 


10  I  Y-1 


T  T  T  -  1 
^20  32  54 


From  Equation  F.9.1 


20 


1  + 

^  ^  2  '^0, 


Solve  for  M^: 


m2  =  M-  T 
”0  \  Y-1  /  20 


Dividing  Equation  F.9.13  by  Equation  F.9.12 


M 

0 


T  T  T  -  1 
20  32  ^54 

T  -  1 
20 


(F 


(B 


(I 


.9.11) 


.9.12) 


.9.13) 


.9.14) 


L 


F.31 


DERIVATI®!  F.9  (Continued) 
Fran  Equation  F.9. 4 


^10  ^10  V^^'^10 


^0  =  ""o 


Hence 


^10  \  ^10 


(F.9. 15) 


where  Equation  F.9. 11  was  used  in  the  last  step. 


Combining  Equations  F.9. 15  and  F.9. 16 

/  \2 


T  T  T 
20  32  54 


\  "^20  “  ^ 


(F.9. 16) 


Since  =  Wn, 
c  T 


^r3  ~  \2  ^4  “  ^5 


T  —  T  =  T  —  T 
T3  T2  T4  T5 


32  ■ 

1) 

=  T, 

1  —  . 

^"^32 

~  1) 

X 

'^43 

^32^ 

(F.9. 17) 


Substitute  Equation  F.9. 17  into  F.9. 16 


(T20  -1)  ^*^43“^^  ^^32 


^'^43^  (F.9. 18) 
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DERIVATICN  F.9  (Continued) 


From  Equation  F.9. 5 


^n  =  ?  ^0  - 


P 

g  0  V, 


(F.9. 19) 


Substituting  Equation  F.9. 18  into  F.9. 19 


^n=|^0  <^43  <'^32  -  +  ^43  -  1 


(F.9. 20) 


Equation  F.9.20  shaTS  that  net  thrust  is  dependent  only  on  three  design 
parameters. 


T.,  =  = —  proportional  to  fuel  flow 

^T3 


T,-  =  = —  proportional  to  coirpressor  ratio 

T2 


proportional  to  Mach 


If  1^2  approaches  1,  a  ramjet  results  and  Equation  F.9.20  reduces  to 


W 

F  =  -  V_ 
n  g  0 


(F.9. 21) 


Iilhat  is  the  static  thrust  (M  =  0)  of  the  ramjet  and  turbojet? 


Ramjet: 


lim  F  =  lim  —  V» 

Mq^O  Mq^O  ^  ^ 
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DERIVATiast  F.9  (Continued) 


Note;  V-  =  a  M„  vdiere  a  is  the  speed  of  sound  and  T.-  is  independent  of  Mach 


This,  of  coiirse,  is  the  ejqjected  result  as  a  ramjet  does  not  produce  any 
static  thrust. 


Turbojet; 


-  .  w 

Lun  F  =  lim  -  V_ 

M  n® 

0 


Tt;^)  '•^43  -1'  «32  -1 


Since 


T  =  T 
T2  TO 


T  =  -12  =  1  +  :LZ-i  m2 

^20  Tq  2  ^0 


T  1  +  ^  ~  ^ 

20  ^  ^  2  % 

^20-1' 


Substituting  and  simplifying 


=  |a  lim 


1  +X^m2  ^ 

Y  -  1  ^*^43  “  ^"^32  "  *^43  ^0  "  ^0 
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DERIVATICX^  F.9  (Continued) 


Thus 


vMch  shows  that  unlike  the  ramjet,  the  tui±»ojet  produces  thrust  at  zero 
velocity. 
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F.IO  ALTERNATE  DERIVATION  OF  IDEAL  NET  THRUST  EQUATION  FOR  A  TURBOJET 
(Chapter  7) 

Starting  with 


i  <^10  -  ''o' 


(F.10.1) 


Wfe  want  to  express  in  terms  of  engine  parameters  and  flight  conditions. 
Fran  cycle  analysis 


=  2gJ  (h^g  -  h^g)  vhere  =  h^g 

=  2gJ  (hpg  -  hg  )  where  h^g  =  h^2 


(F.10.2) 


(F.IO. 3) 


Siobtracting  Equation  F.IO. 3  from  Equation  F.10.2 


2  2 
^10  “  ^0  = 


2gJ  1^5  -  hj^g  -  11^2 


(F.IO. 4) 


Since 


-  lip3  -  lip2  '^^T  ^4  "  ^5 


“  ^2  ^4  “  ^3 


(F.IO. 5) 


Substituting  Equation  F.10.5  into  Equation  F.IO. 4 


■  ''o  “  ^  [^4  ”  ''lo'  '  **'r3  ■  *'o'j 

=  2gJCp[iT^^  -  Ti„)  -  (Tj3  -  T^)] 


=  "t4 


(F.IO. 6) 
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DERIVATIOJ  F.IO  (Continued) 


However  for  isentropic  flow 


Y  -  1 


(F.10.7) 


"^0  V^IO 


Y  -  1 

Y 


Y  -  1 


Tt4  J^tA  ^  f^T3\  ^ 


(F.10.8) 


since  for  an  ideal  process 


P  =  P 
T3  ^T4 


'■lO  =  •’o 


Substituting  Equations  F.10.8  and  F.10.7  into  Equation  F.IO. 6  and  subtracting 

r  _  ^ 


Y  -  1  Y  -  1 

J  /P  \  ^  /P  \  ^ 

=  2gJC  1  -  (—  1  -  T  (  -^1 

10  ^  p  j  T4  j 


-  1  rv: 


(F.10.9) 


P  P 

T3  ^  ^T3 

P  P 

0  ^T2 


1 

^T1 

•  •  :^  =  CR  (1  + 

7T1  ^0 


Y  -  1  w  . 

Mrt) 

Y  0 


=  CR  f  (M) 


The  last  term  is  the  ram  recovery,  f  (M) 
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DERIVATION  F.IO  (Continued) 

Solving  for  =  TIT  in  Equation  F.10.9, 


=  f  (TIT,  CR,  Mq,  Tq) 


F.38 


